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The new observational data on interstellar dust, obtained during the past few years, in particular, revised 
ISM elemental abundances, the polarization curve from the far UV to the near IR, spectra of dust infrared 
emission, extended red emission (ERE), are in serious conflict with existing dust models. We proposed 
here to use the regularization approach as a base for developing a modern dust model, which attempts to 
self-consistently explain most of the available observational data. 



1. Why is dust important for astronomy? 

Dust plays an important role in various cosmic environments: circumstellar shells, star-forming 
regions, protoplanetary nebulae, interstellar and intergalactic medium (for a review see Dorschner 
k Henning |l99^ ). The dust grains absorb, scatter and re-emit the electromagnetic radiation from 



various cosmic objects. The dust grains contribute to the thermal equilibrium of interstellar 
matter. Various chemical and physical processes, e.g. formation of molecules, take place on the 
surfaces of dust grains. Thus, the knowledge of the properties of cosmic dust is important for: 
1) the correct interpretation of the observed spectra of objects seen through intervening dust 
and the observations of radiation emitted from, or scattered by dust and 2) modeling of objects 
where dust is a major constituent. 

We may mention a few dust models which were most popular until recently. One of them 
is the so-called MRN model (Mathis, Rumple & Nordsieck 1977] ): a mixture of spherical bare 



graphite and silicate grains, having a power law size distribution /(a) ~ a over the range 
of sizes 0.005 - 0.25 /im. This model was further developed by Draine & Lee ( [1984| ), Mathis &: 
Whiffen (|1989|) , Kim, Martin & Hendry (|1994| ), Mathis (|1996| ). Another dust model was proposed 



by Hong & Greenberg ( |1978D : small bare graphite grains (for explanation of the 217.5 nm hump) 
and silicate grains covered by the organic refractory mantles produced by UV photolysis of dirty 
ices which are formed on the grains in molecular clouds. It was assumed that such complicated 
silicate grains obey some exponential cubic size distribution. This model was further developed 
by Greenberg (|l989| ), Greenberg & Li ( [1996D , Li & Greenberg ( |l99^ ). Both ab ove-mentioned 



models were in rather good agreementwith the mean Galactic extinction and polarization laws, 
provided the cosmic elemental abundances are equal to solar abundances. 

2. New observational data on cosmic dust 

During the past few years considerable progress has been achieved in obtaining new observational 
data on interstellar dust. Here are some of the most important findings: 

(a) It became clear that the chemical compositions of the general interstellar medium and 
the Sun are not the same, as was previously assumed: the cosmic elemental abundances 
appear to be 60-70 per cent of the solar abundances (Snow & Witt 1996|) . 



(6) The analysis of high signal-to- noise spectra of interstellar carbon and oxygen in a few 
directions obtained with the Goddard High Resolution Spectrograph aboard the Hubble 
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Space Telescope, produced the first reliable estimates of the mean gas-phase abundances of 
interstellar carbon: 10^C/HRil40±20 and oxygen 10^O/H?a310±20 (Cardelh et al. |1996D . 



(c) It was found that the dust albedo in the near-IR, 0.6-0.8, considerably exceeds the 
albedo predicted by the standard dust models (Witt et al. 1994 ; Lehtinen & Mattila |1996|) . 
Recently, Witt et al. ( 1997| ) derived estimates of both albedo and asymmetry parameter of 
dust in the diffuse medium in the far-UV, which are generally consistent with expectations. 



(d) Far-UV polarization observations have been performed (Clayton et al. 1992; Somerville 
et al. 1994 ; Clayton et al. 1995| ), and produced the wavelength-dependent interstellar po- 



larization curves from the far-UV down to the near-IR. 

(e) The infrared emission from the diffuse interstellar medium in the range 3.5-1000 fim 



has been measured for the first time with COi?-E'/DIRBE and FIRAS (Bernard et al. 1994; 
Dwek et al. |1997| ; Sodroski et al. |1997l ; Lagache et al. |1998|) . 



(/) The so-called extended red emission (ERE) in the 500-800 nm spectral range has been 



detected in the general interstellar medium (Gordon et al. 1998; Szomoru & Guhathakurta 
|1998| ). Witt et al. (|1998| ) and Ledoux et al. (1998) independently proposed that silicon 
nanoparticles might be the source of the ERE. 



3. Previous work on modeling of interstellar dust 



With the wealth of new data it became clear that practically all dust models, proposed before 
1995 (Dorschner Sz Henning |1995| ), face a serious crisis in explaining the new observational data. 
No attempts were undertaken so far to build a dust model which would accomodate most of 
the data recently acquired. However, some important theoretical steps have been done towards 
such the model. In particular, Mathis ( 1996| ), Zubko et al. ( 1996| , 1998) and Li & Greenberg 
( 1997 ) proposed more sophisticated dust models on the base of composite, core-mantle and 
multilayer grains with the main goal to explain the interstellar extinction using the revised 



elemental abundances. Li & Greenberg ( 1997 ) have made an attempt to analyse the interstellar 
extinction and polarization, taking into account the new chemical abundances. Dwek ( 1997] ) has 
tried to model the extinction and infrared emission together. Note that the approaches used by 
Mathis (HH), Dwek ([T997| ) and Li & Greenberg (|l997| ) were based on the predefined grain-size 
distributions: power law or Gaussian and, therefore, cannot be considered as leading to the best 
solution. Kim et al. (1994) and Kim & Martin (|1994 |1995| , |1996D have modelled the interstellar 
extinction and polarization using the Maximum Enthropy Method (MEM). However, the authors 
treated separately the extinction and the polarization rather than self-consistently. In addition, 
as was noted by Zubko et al. ( 1996| ), the MEM cannot produce a unique solution, because some 
default solution should be specified prior to the calculations. As a rule, the default is not known 
a priori and this introduces some uncertainty into the final result. 

Recently, Zubko ( 1997] ) has proposed the regularization approach for modeling of interstellar 
dust. To the best of our knowledge, it is the only tool which is capable of deriving an optimum 
and unique grain-size distribution in a general form for any predefined mixture of model grains. 
In addition, the regularization approach directly takes into account the mathematical nature of 
the problem, being an ill-posed inverse problem. A few years ago the author started to develop 
a computer program based on the regularization approach. In its present state, this program is 
capable of building a model of interstellar dust by a self-consistent analysis of an extinction curve, 
scattering properties: albedo and asymmetry parameter, elemental abundancies and mass fraction 
constraints. Currently, the program works with spherical grains which may be homogeneous (their 
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optical properties are calculated with the standard Mie approach), multilayer (author's original 
approach) or of composite structure (EMT/Mie approach). The first models of both circumstellar 
and interstellar dust derived with the program appear to be very promising (Zubko [1997| ; Zubko 



et al. 199e, 1998). 



4. What should a modern dust model fulfill? 

In our view, an ideal model of interstellar dust, which could be constructed in near future, should 
fulfill the following requirements: 

(a) it should be a unified dust model, that is, it should accomodate most existing observa- 
tional data on extinction, polarization, elemental abundances and so on. 

(6) it should use dust particles of more realistic structure: composites, fractal, core-mantle 
and multilayer grains, both spherical and non-spherical; 

(c) it should allow more grain species, and use the optical properties, preferably based on 
laboratory measurements of cosmic analogue grains; 

(d) it should handle the grain-size distributions in the most general form rather than in 
some prescribed form, e.g. a power law or a Gaussian; this approach should provide an 
optimum solution. 

The best-fit model will be derived using the regularization approach, which in addition to the 
already mentioned advantages, will allow one to easily include other observational constraints in 
the future. One of the by-products of this work will be a well-documented program for generating 
a self-consistent model of interstellar dust, which will be available upon request. 
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